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ABSTRACT 

We study the detection error probability associated with bal- 
anced binary relay trees, in which sensor nodes fail with 
some probability. We consider N identical and independent 
crummy sensors, represented by leaf nodes of the tree. The 
root of the tree represents the fusion center, which makes the 
final decision between two hypotheses. Every other node is a 
relay node, which fuses at most two binary messages into one 
binary message and forwards the new message to its parent 
node. We derive tight upper and lower bounds for the total 
error probability at the fusion center as functions of N and 
characterize how fast the total error probability converges 
to with respect to N. We show that the convergence of 
the total error probability is sub-linear, with the same decay 
exponent as that in a balanced binary relay tree without sen- 
sor failures. We also show that the total error probability 
converges to 0, even if the individual sensors have total error 
probabilities that converge to 1 /2 and the failure probabilities 
that converge to 1, provided that the convergence rates are 
sufficiently slow. 

Index Terms — Binary relay tree, crummy sensors, dis- 
tributed detection, decentralized detection, hypothesis testing, 
information fusion, dynamic system, invariant region, error 
probability, decay rate, sensor network. 

1. INTRODUCTION 

Consider the decentralized detection problem introduced 
in 0]: Each sensor makes a measurement and summarizes its 
measurement into a message. These messages are forwarded 
to the fusion center, which then makes a final decision. 

This decentralized detection problem has been studied in 
the context of several different network topologies. In the 
parallel architecture, also known as the star architecture [T]- 
Ifl5ll . ll3"2l . all sensors directly communicate with the fusion 
center. When sensor measurements are conditionally inde- 
pendent, the decay rate of the total error probability in the 
parallel architecture is exponential Q. 

This work was supported in part by AFOSR under Contract FA9550- 
09-1-0518, and by NSF under Grants ECCS-0700559, CCF-0916314, and 
CCF-1018472. 



Another well-studied configuration is the tandem network 
1 1 6 1 — 1 20 1,| 32]. The decay rate of the error probability in this 
case is sub-exponential GUI . Furthermore, as the number of 
sensors N goes large, the error probability is Vl(e~ cN ) for 
some positive constant c and for all d £ (1/2, 1) [18 |. This 
configuration represents a situation where the length of the 
network is the longest possible among all networks with N 
nodes. 

The configuration of bounded-height tree has been studied 
in |21 1-|29|,|32|. This configuration reduces the transmission 
cost compared to the parallel configuration. In the bounded- 
height tree structure, leaf sensor nodes summarize their mea- 
surements and send the new messages to their parent nodes, 
each of which fuses all the messages it receives with its own 
measurement (if any) and then forwards the new message to 
its parent node at the next level. This process takes place 
throughout the tree culminating in the fusion center, where 
a final decision is made. If only the leaf nodes are sensors 
making measurements, and all other nodes simply fuse the 
messages received and forward the new messages to their par- 
ents, this tree is known as a relay tree. For a bounded-height 
tree with lim rjv _ i . 00 £n/tn = 1, where tn denotes the total 
number of nodes and denotes the number of leaf nodes, 
the optimum error exponent is the same as that of the parallel 
configuration J22). 

For trees with unbounded height, the convergence anal- 
ysis is still largely unexplored. In [30], the convergence of 
the total error probability in balanced binary relay trees with 
unbounded height has been proved. Upper and lower bounds 
for the total error probability at the fusion center as functions 
of TV have been derived in [31 1. These bounds reveal that the 
convergence of the total error probability at the fusion center 
is sub-linear with a decay exponent ^/~N. 

In this paper, we assume that each of the sensors fails with 
a certain probability. A failed sensor will not provide a mes- 
sage to its parent node at the next level. We refer to these 
sensors as crummj^ sensors. We will derive upper and lower 
bounds for the total error probability at the fusion center as 
functions of N. Not surprisingly, we find that the decay of 

'The attentive reader will recognize that our use of the term 
"crummy"follows in the footsteps of our great patriarch, Claude E. Shannon 



the total error probability for each step is worse than the case 
where there is no sensor failure. But this decay rate is still 
sub-linear with the same decay exponent vN in the asymp- 
totic regime, regardless of the sensor failure probability. 

2. PROBLEM FORMULATION 

We consider the problem of binary hypothesis testing between 
Ho and Hi in a balanced binary relay tree with crummy sen- 
sors. As shown in Fig. [T] leaf nodes are sensors undertak- 
ing initial and independent detections of the same event in a 
scene. These measurements are summarized into binary mes- 
sages. If a sensor node works properly, then it forwards the 
summarized message to its parent node at the next level. Oth- 
erwise, with a certain probability the sensor fails in the sense 
that it does not forward the message upward. Each non-leaf 
node — except the root, which is the fusion center — is a relay 
node, which fuses binary messages it receives (if any, and at 
most two) into one new binary message and forwards the new 
binary message to its parent node. This process takes place 
at each intermediate node culminating the fusion center, at 
which the final decision is made based on the information re- 
ceived. 




that the two child nodes have identical Type I error probability 
a, identical Type II error probability /3, and identical failure 
probability q. 

Denote the Type I error, Type II error, and failure proba- 
bilities after the fusion by (a' , f3' ,q'). This parent node fails 
to provide any message to the node at the next level if and only 
if both its child nodes fail to forward any message. Hence, we 
have 

q' = q 2 - (1) 

If one of the child nodes fails and the other one sends 
its message to the parent node, then Type I and Type II er- 
ror probabilities do not change since the parent node receives 
only one binary message. The probability of this event is 
2q(l — q), in which case we have 



(2) 



If both child nodes send their messages to the parent node, 
then the scenario is the same as that in |30| and [31]. The 
probability of this event is (1 — q) 2 , in which case we have 



(l-(l-a) 2 ,/3 2 ), a</3, 
(a 2 ,l-(l-/3) 2 ), a>p. 



(3) 



Let a' and j3' be the mean Type I and Type II error proba- 
bilities conditioned on the event that at least one of these child 
nodes forwards its message to the parent node, i.e., the parent 
node has data. We have 



(a 2 ,/3 2 ,q 2 ) 

(«i.A>9i) 
(a ,p ,q ) 





(a',P',q') = f(a,p,q) 



(4) 



N 



Fig. 1. A balanced binary relay tree with height k. Circles 
represent sensors making measurements. Diamonds represent 
relay nodes which fuse binary messages. The rectangle at the 
root represents the fusion center making an overall decision. 

We assume that all sensors are independent given each hy- 
pothesis, and that all sensors have identical Type I error prob- 
ability o;o and identical Type II error probability (3q. More- 
over, we assume that all sensors have identical failure prob- 
ability q Q . Assuming equal prior probabilities, we use the 
likelihood-ratio test 1 33 ] when fusing binary messages at in- 
termediate relay nodes and the fusion center. 

Consider the simple problem of fusing binary messages 
passed to a node by its two immediate child nodes. Assume 



(l-q)(2a-a 2 )+2qa (l-q)P 2 +2q/3 n i 
1+9 ' 1+9 



q 2 ),a<(3, 

( 

,<7 2 ) ,a > fi. 



(l-q) a 2 +2qa (l-q)(20-l3 2 )+2q0 ^ 2 
1+9 ' 1+9 



Our assumption is that all sensors have the same er- 
ror probabilities (ao, (3q, <7o)- Therefore by |5]), all relay 
nodes at level 1 will have the same error probability triplet 
(«!,/?!, qx) — /(ao,/?o,9o) (where a\ and fli are the con- 
ditional mean error probabilities). Similarly by (4), we can 
calculate error probability triplets for nodes at all other levels. 
We have 



(afe+i^fe+i^fe+i) = f(a k ,f3 k ,q k ), k = 1,2, 



(6) 



where (a k ,Pk, Qk) is the error probability triplet of nodes at 
the fcth level of the tree. Notice that if we let q — 0, then the 
recursive relation reduces to the recursion in OTI . 

The relation ^ allows us to consider (a k , Pk,qk) as a 
discrete dynamic system. For the case where qo = 0, we 
have studied (See Olll ) the precise evolution of the sequence 
{(&k, Pk)}, derived total error probability bounds as func- 
tions of N, and established asymptotic decay rates. In this 



paper, we will study the case where qo 7^ 0. We will derive 
total error probability bounds and determine the decay rate of 
the total error probability. 

To develop intuition, let us start by looking at the sin- 
gle trajectory shown in Fig. |2ja), starting at the initial state 
(ao, Po, qo)- We observe that q k decreases very fast to 0. In 
addition, as shown in Fig. [2jb), the trajectory approaches 
/3 = a at the beginning. After (ak,f3k) gets too close to 
j3 = a, the next pair (ctfe+i, Pk+i) will be repelled toward 
the other side of the line ft = a. This behavior is similar to 
the scenario where q = 0. For the case where q — 0, there 
exist an invariant region in the sense that the system stays in 
the invariant region once the system enters it [31 1. Is there an 
invariant region for the case where q ^ 0? We answer this 
question by precisely describing this invariant region in R 3 . 




a 



(b) 

Fig. 2. (a) A typical trajectory of (a k , (3 k ,q k ) in the (a, /?, q) 
coordinates, (b) The trajectory in (a) projected onto the (a, j3) 
plane. 



3. THE EVOLUTION OF TYPE I, TYPE II, AND 
SENSOR FAILURE ERROR PROBABILITIES 

The relation |5]l is symmetric about the hyper-planes a + ft = 
1 and f3 — a. Thus, it suffices to study the evolution of the 
dynamic system only in the region bounded by a + f3 < 1, 



P > a,and0 < q < 1. LetU := {(a, (3) > 0|a+/3 <l,(3> 
a, and < q < 1} be this triangular prism. Similarly, define 
the complementary triangular prism C := {(a,/3) > 0\a + 
P < 1,jS < a, and < q < 1}. 

First, we denote the following region by Si :— { (a, (3, q) e 
U\& < (-q+^q 2 + (1 - q) 2 (2a - a 2 ) + 2q(l - q)a)/(l- 
q)}.\f(a k ,l3 k ,q k ) € B x , then the nextpair (a k +i, (3 k +i, q k +i) 
jumps across the plane j3 = a away from (a k , f3 k ,q k ). More 
precisely, if (a k ,f3 k ,q k ) E U, then (a k ,/3 k ,q k ) £ B x if and 
only if (a k+ i, f3 k+ i,q k+ i) <E C. This set Bi is identified in 
Fig.^a). 

It is easy to see from <|3j and (|6]) that, if we start with 
(ao, /3o, qo) G U \ B\, then before the system enters B\, we 
have Qffc+i > a k and P k +\ < (3 k . Thus, the system moves 
towards the j3 — a plane. Therefore, if the sensor number 
is sufficiently large, then the system is guaranteed to enter 
B x . 




Fig. 3. (a) Region B\ in the (a, (3, q) coordinates, (b) Region 
Ru in the (a, (3, q) coordinates. 

Next we consider the behavior of the system after it enters 
B\. If (a k , (3 k ,q k ) € B\, we consider the position of the 
next pair (afc+i, f3 k +i, q k +i), i.e., consider the image of Bi 
under /, denoted by i?£. Similarly we denote the reflection 
of Rc with respect to /3 = a by Ru. We find that Ru :— 



{(a,/3,q)€U\p < -a + 2(V? + (1 - q 2 )a- q)/(l- q)} 
(seeFig.|3Jb)). 

The sets i?^ an d Si have some interesting properties. We 
denote the projection of the upper boundary of Ru and B\ 
onto the (a, f3) plane for a fixed q by R u and B^, respectively. 
It is easy to see that if qi < q2, then Rff lies above R^f in 
the (a,/3) plane. Similarly, if qi < q 2 , then Bf 1 lies above 
Bf 2 in the (a, /3) plane. Moreover, we have the following 
Proposition. 

Proposition 1: B\ C i?^. 

Proof. B\ and i?^ share the same lower boundary f3 = a. 
Thus, it suffices to proof that the upper boundary of B\ is 
below that of Ru for a fixed q, i.e., R u lies above B^ in the 
(a, /?) plane. 

The upper boundary of B\ is 

p = -g + yV + (i - g) 2 (2« - « 2 ) + 2g(i 

1-9 

The upper boundary of Ru is 

^ = - a + 2^ 2 + (1 '' 72)a ^. 

i - q 

Notice that when q — 0, these boundaries reduce to the 
boundaries in JST). We need to prove the following: 

-q + yV + (1 - g) 2 (2a - a 2 ) + 2g(l - q)a 

i_-q 

< - Q + 2^ 2 + (1 - g2)Q ^. 

1 - 5 

It suffices to show that 

vV + (1 - g) 2 (2a - a 2 ) + 2g(l - q)a 

< -a(l - g) - 9 + 2vV + (1 - q 2 )a. 

Squaring both sides and simplifying, we have 

2y/q 2 + (l-q 2 )a(a(l - q) + q) 

< 2(q 2 + (l-q 2 )a)-(l-q) 2 (a-a 2 ). 

Again squaring both sides and simplifying, we have 

4(q 2 + (1 - q 2 )a){q 2 + 2q(l - q)a + (1 - qfa 2 
-q 2 -(l-q 2 )a + (l-q) 2 (a-a 2 )) 
< {l-q)\a-a 2 ) 2 . 

Fortuitously, the left hand side turns out to be identically 
0. Thus, the inequality holds. The reader can refer to Fig. |4|a) 
and Fig. |4|b) for plots of the upper boundaries of Ru and B\ 
for two fixed values of q. 

□ 
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Fig. 4. (a) Upper boundaries for Ru and B\ for q = 0.1. (b) 
Upper boundaries for Ru and £?i for q = 0.01. 

We denote the region Ru U i?£ by R. We show below 
that R is an invariant region in the sense that once the system 
enters R, it stays there. 

Proposition 2: If (ot-koi Pk a ,<lk ) £ i? for some fco, then 
{a k ,l3 k ,q k ) £ R for all > k . 

Proof. Without lost of generality, we assume (a k , /3fe, qk) G 
We know that i?£ is the image of U in £. Thus if the 
next state (a k +i, Pk+i, 9fe+i) G then it must be inside 
We already have q k +i < <7fc, which indicates that i?^. fc+1 
lies above R^ in the (a, (3) plane. Moreover, for a fixed q, 
the upper boundary R u is monotone increasing in the (a, f3) 
plane. We already know that a k +i > a k and fik+i < Pk- 
As a result, if the next state (a k +i, Pk+i, Qk+i) G U, then the 
next state is in fact inside R u . 

□ 

We have shown that the system enters Bi after certain 
levels of fusion. By the fact that B\ C Ru, we conclude that 
the system enters Ru at some level of the tree and stays inside 
the invariant region R at all levels above. 

In the next section, we will consider the step-wise reduc- 
tion of the total error probability when the system lies inside 



the invariant region and deduce upper and lower bounds for 
the total error probability. 

4. ERROR PROBABILITY BOUNDS 

The total detection error probability for a node at the kth level 
is (a k + Pk) l '2 because of the equal -prior assumption. Let 
Lk = a k + Pk, which is twice the total error probability. We 
will derive bounds on logL^" 1 , whose growth rate is related 
to the rate of converge of L k to 0. (Throughout this paper, log 
stands for the binary logarithm.) 

Proposition 3: Let L^}^ be the total error probability at 
the next level from the current state {ot k ,Pk,q)- Suppose that 
(a k , Pk,qi) and (a*., P k , q 2 ) £ U. If qi < q 2 , then 

L ( , q l] < L (q2) 



J k+1 



Proof. From we have 



-(<?) 

J k+1 



l + « 



q L° 



fc+r 



2q_ 
l+q 



(a k + Pk), 



where = 2a k - oi{ t Hk - 

It is easy to show that 2a k — a\ + Pi < a k + Pk. 

2a k -a 2 k + pl<a k + p k 
a k - a\ < f3 k - Pi. 

Since a k + p k < 1 and p k > a k , we have p k — 1/2 < 
1/2 — a k . Notice that the function x — x 2 peaks at x = 1/2. 
Hence, 2a k - a\ + Pi < a k + p k . 
Notice that 

1-q , 2q 



Pk 



l + q l + q 
Therefore we can write 



= 1. 



L k %{=PiL° k+1 + (l- Pl )(a k +p k ), 

where pi = (1 - qi)/(l + qi). Let p 2 = (1 - £ya)/(l + fe), 
it is easy to see that pi > p2- Thus we have 



L 



in) 

k+l 



= PiL° k+1 + (1 - Pl )(a k + p k ) 

+(P2 - Pi)L° k+1 - (p 2 - Pi)L° k+1 
< Pl L° k+1 + {1 - Pl )(a k + p k ) 

+(P2 -pi)L° k+1 - (p 2 - P i)(a k + P k ) 



= L 



(92) 
fe+i- 



□ 



From Proposition 3, we immediately deduce that 



L 



(o) 
k+l 



< L 



(91) 
fe+i- 



This means that the decay of the total error probability for a 
single step is the fastest when q = 0. As a result, for the 



case where g / 0, the step-wise shrinkage of the total error 
probability cannot be faster than the case where q — 0, where 
the asymptotic decay exponent is \/N OTI . 

Notice that from ([TJ, the decay of q k is quadratic, which is 
much faster than the decay rate of L k . Moreover, it is easy to 
see that the decay of q k is faster than the decay of a k and of 
P k . Hence, it is natural to assume that q k < a k and q k < p k 
when we consider the step-wise shrinkage of the total error 
probability in the invariant region. Next we give upper and 
lower bounds for the ratio L k+2 /Ll. 

Proposition 4: Suppose that (a k , Pk,qk) & R, OL k > q k , 
and p k > q k . Then, 

Proof. First, we consider the lower bound. The evolution of 
the system is 

(a k ,p k ,q k ) -> (a k+1 ,p k+1 ,ql) -> (a k+2 , (3 k+2 ,q k ). 
From Proposition 3, we have 

r(0) < j 
J^ k+2 S ^k+2, 

where Lj.°? 2 = ^ a k+i — a k +i + Pk+i as defined before. 
To prove 1/2 < L k + 2 /L^, it suffices to show that 1/2 < 

If (a k ,p k ) e i? u \B 1 ,then 



L 



(o) 

fc+2 



'k+l 



(a k + p k y 



We have 

ctk+i = 



1 - qk rn _ „ 2 N , 2 ?fe 



l + qk 



{2a k - at) + 



l + qk 



ak > a k 



and 



2g fc 

l + qk rK l + qk 
Thus, it suffices to show that 



2a k ~al+Pi 1 
{a k +p k ) 2 -2' 

It is easy to see that 

2{2a k -a 2 )>l-{l-a k )\ 
Hence, it suffices to show that 

(l-(l-a k ) 4 + p 4 k )>(a k +p k ) 2 , 



which has been proved in [31]. 

If (a k , P k ) £ B\, then it suffices to show that 



*fe+l + 2Pk+l - P k +i 1 



{a k + p k f 



We have 



Hence, we have 



"fe+i = (2a k - a fe J + — a k > a k 

1 + Qk L + Qk 



and 



« - 1 ~ qk ft 2 l 2<lk R > « 2 
1 + Qk t + Qk 



Thus, it suffices to proof 



> 



K + /3 fc ) 2 " 2 

which is easy to see. 
Next we prove th 
If (a fe ,/3 fe ) e Si, then 



Next we prove the upper bound of the ratio Lk+ijL\. 



Lk+2 < Lk+\ 



l-q k 2a k -a 2 + fi 2 2q, 



l + qk (afc + Pfc) 2 l + (?fe 
It is easy to see that 

^-{OL k +P k )<l. 

1 + Qk 
Next, we can prove that 

K + / 3 fe )2 " ' 

which is equivalent to 

4>{a k ,p k ) := 2a k - 3a 2 -ft- 4a k L3 k < 0. 
We have 



dcf) 
Wk 



-2(3 k - 4a k < 0. 



Thus, we can consider the lower boundary of this region 
which is the upper boundary of B\. 



-Q + yV + (1 ~ q) 2 (2a - a 2 ) + 2g(l - g)a 
l-q 



Denote tp(a,q) := ^q 2 + (1 - <?) 2 (2a - a 2 ) + 2q(l - q)a. 
We have 

<t>(a k ,P k ) = -(q 2 k + q 2 k + (l-q k ) 2 (2a k -a 2 k ) 

+2q k (l - q k )a k - 2q k ip(a k , q k ))/ (1 - q k ) 2 
-Aa k f3 k + 2a k - 3a 2 . 
2q k Pk . a 2q k a k 2 



1 - 9fe 
It is easy to see that 

2q k /3 k 
1 - qu 



1 - qk 



4a k p k < 0. 



and 



1 - g k 2a k - a 2 + 0j < ^ 
l + 9fc (afc+Pfc) 2 



< 3. 



For the case where (o^, e Si, we prove that the ratio 
is upper bounded by 4. The evolution of the system is 

(a k ,/3 k ,q k ) -> (afc+i^fe+i,^) -> (a k+2 , /3 k+ 2,qt)- 
It is easy to see that 

where 4+ 2 denotes the total error probability if we use qk to 
calculate from L k+ i to L k+2 . Therefore, it suffices to prove 
that 

4+2 - 4i fc = «fe+2 + h+2 - 4(a fe + h) 2 < 0. 
We have 

h + , = \-^Pl + ^Pk. 
l + q k l + q k 

From the assumption that j3 k > q, we have 

9^ = 2(1^) _2«fc_ 
<9pfc 1 + g fe 1 + ft 

It is easy to get that 

Pfc+2 = gfc (2/3 fc +i — /9fc + i) + -i qk Pfc+i 
J- + ft J- + ft 



1 - ft o2 

l + ft Pfc+1 ' l + qk 



'-0t+i + ——Pk+i- 



Therefore, we have 



1 - Qk R d/3 k +i 2 dh+i Qfl 

op k l + qk o(3 k 1 + q k d(3 k 

Thus, 

dL[ qk i - AL 2 

Therefore, we can consider the lower boundary of B\, j3 k — 
a k . We have 

r (m) , f 2 _ 4(l-g fc ) 2 (l-g fc ) 2 , (l-gfc) 2 3 



+ 



2(1 -ft) 2 2 



(1 + ft) 2 " fc+ (l + q k ) 2 



a k - 16a fe < 0, 



which holds in region Si. Hence, the ratio is upper bounded 
by 4 in this region. 

□ 



Proposition 4 gives rise to bounds on the change in the 
total error probability every two steps: L k+2 < and 
Lk+2 > L\/2. From these, we can derive bounds for log LIT 1 
for even-height trees, i.e., fc = log TV is even. Let Pjy — 
L\ og N, namely, the total error probability at the fusion center. 
We will derive bounds for log P^ 1 . 

Theorem 1. If(a , f3 , q Q ) e R and log TV is even, then 

Proof. If (a ,(3 ,q ) G R, then we have (a k ,/3 k ,q k ) € R 
for k = 0, 1, ... , log TV — 2. From Proposition 4, we have 

for A: = 0, 1, ... , log TV — 2 and some e [1/2, 4]. There- 
fore, for k = 2, 4, . . . , log TV, we have 



r fe/2 

n 



r.k/2 



where a,- L e [1/2, 4]. Substituting fc = log TV, we have 



'log v^V 



p N =\ n °* 



r2'° B ' 



i=l 




Hence, 



/7v 



logP^ = - ^ log 



i=l 



Notice that logi 1 > and for each i, — 1 < loga 4 < 2. 
Thus, 

logP^ 1 < y/NlogL^ 1 +log\/iV 

log v 7 ^ \ 



Finally, 



TV log L^ 1 



-l 



log P^ 1 > -21og VTV + VN log L 

21og ViV 



VN logL^ 1 



□ 



For odd-height trees, we need to calculate the decrease in 
the total error probability in a single step. For this, we have 
the following Proposition. 



Proposition 5: If (a k ,(3 k , qk) ^U, then we have 
Lk+i 



and 



T 2 
L k 



Lk 



> 1 



< 1. 



Proof. To prove L fc+1 /L| > 1, it suffices to prove that 



1 - qk 



l + qk 

+ 2qk 



(2a k -a 2 k + pi- K + p k ?) 



l + qk 



{a k +fi k - (a k +/3 fe ) 2 ) > 0, 



which is easy to see. 

To prove L k+ i/L k < 1, it suffices to prove that 



i - qk 



1 + qk 

+ 2q k 



{2a k ~al+fil-{a k +p k )) 



i + qk 

which is easy to see. 



K+At- (a k + p k )) < 0, 



□ 



From Propositions 4 and 5, we give bounds for the total 
error probability at the fusion center for trees with odd height. 



Theorem 2. If(a , /3 , <7o) € R> then 



log^ 1 - 



2 log < > * 



< %/2TV I logig 1 



Proq/\ By Proposition 5, we have 



^r2 



< - log Pat 



logyf 

\/2TV 



ii = aL 



for some a > 1. 

By Proposition 4, we have 

= 



a k L\ 



for fc = 1,3, .. . , log TV — 2 and some a k e [1/2,4]. Hence, 
we can write 



'(k-l)/2 

ife = a | ] | ai 



-,(fc + l)/2 



where 1/2 < a, < 4 for i = 1, 2, . . . , (fc - l)/2 and a > 1. 
Let fc = log TV, we have 



P N = a\ H aA^™ =~a\ H aALf*, 



log 



i=l 



and so 



log -P^ 1 = -logo- 



log en I +V2NlogL^\ 

V 4=1 



Notice that logL > and for each i, log a,i > —1. More- 
over, log a > 0. Hence, 

TV 



logP^ 1 < V^logL^+logW- 



= V2N^logL 

By Proposition 5, we can write 

Li = aLo 

for some a < 1. Thus, 



io g yf 

i/2TV 



'(fc-l)/2 

Lfe = a | J~| 



2 (fc-l)/2 



where 1/2 < a f < 4 for « = 1, 2 . . . , (Jfe - l)/2 and a < 1. 
Hence, 



Pjy = a 



n 



V 



V 



and so 

log-P^ 1 = -log 5 - I l°g a * 



n ai i l o 



login? 1 . 



Notice that log L 1 > and for each i, —1 < log < 2 
and log a < 0. Thus, 

logP^ 1 > -21og^- + ^-logL - 1 



log^c 



2 log Jf 



5.1. Asymptotic decay rate 

Notice that when TV is very large, the sequence (3t,qk)} 
enters the invariant region R at some level and stays inside 
afterward. Therefore the decay rate in the invariant region 
determines the asymptotic rate. Because our error probability 
bounds for odd-height trees differ from those of even-height 
trees by a constant term, without lost of generality, we will 
consider trees with even height to calculate the decay rate. 
Proposition 6: If Lq = «o + A is fixed, then 



logP^-logLoW. 

Proof. If Lq — cto + A is fixed, then by Theorem 1 we im- 
mediately see that P/v 
and 



1 - 2 



< 



logP 



A' 



logVTV 
logLj W ~ logL-W 



as TV -> oo (log Pjy -> oo) 
logVTV 



< 1 



login" ViV 



In addition, because log VTV/v TV — ¥ 0, we have 
logP^ 1 k , 



which means 



logL^W 



logP^ 1 ~logLo W. 



□ 



This implies that the convergence of the total error prob- 
ability is sub-exponential with decay exponent \fN. Com- 
pared to the decay exponent for the case where q = (no 
sensor failures), the asymptotic rate does not change when 
we have crummy sensors, even though the step-wise shrink- 
age for the crummy sensor case is worse. 

Given L £ (0, 1) and e £ (0, 1), suppose that we wish to 
determine how many sensors we need to have so that Pm < 
e. The solution is simply to find an TV (e.g., the smallest) 
satisfying the inequality 



> — loge. 



N login 



2 log VN 



TV 



□ 



The smallest TV grows like 0((log e) 2 ) (cf., |5T| , in which the 
growth rate is the same, and BUI , where a looser bound was 
derived). 



5. ASYMPTOTIC RATES 

In this section, we first consider the asymptotic decay rate of 
the total error probability with respect to TV. We compare the 
rate with that of balanced binary relay trees without sensor 
failures. Then we allow the sensors to be asymptomatic ally 
bad, in the sense that qo — > 1 and ao + f3 — > 1. We prove 
that the total error probability still converges to provided the 
convergence of qo and ao + A is sufficiently slow. 



5.2. Asymptotically bad sensors 

First we consider the case where q depends on TV (denoted 
by ?q )• We wish to have the failure error probability at the 
fusion center qi og jy to converge to 0. 

If q^ is bounded by some constant q £ (0, 1) for all TV, 

then clearly q\ og N — > 0. So henceforth suppose that — > 
1, which means that the sensors are asymptotically arbitrarily 
unreliable. 



Proposition 7: Suppose that q = 1 — yyjv with tjn — > 0. 
Then, gi og at — > if and only if i]n = u>(l/N) (i.e., 77^^ — > 
00). 

Proof. From ([TJ, we have 

Ik — Wo J ' 

Letting /c = log TV, we can write 



or equivalently, 

It is easy to see that q\ og N — > if and only if iVlog(l — 
?7jv) _1 — > 00. But as x — > 0, — log(l — x) ~ x/ln(2). 
Hence, giogw — > if and only if t^jv^V — > 00. 

□ 

Now suppose that ci/7V < r/N < C2/N. In this case, for 
large N we deduce that 

Cl < log^gw < C 2 , 

or equivalently, 

2~ C2 <g log iv<2- Cl . 

Finally, if -q^ = o(l/N) (i.e., 777V converges to strictly 
faster than 1/N), then gi og jv — > L 

Next we allow the detection error probability of individual 
sensors to depend on N, denoted by Lq . 

If Lq N ^ is bounded by some constant L g (0, 1) for all 
N, then clearly Pjy — > 0. It is more interesting to consider 
Lq 1, which means that sensors are asymptotically bad. 

Proposition 8: Suppose that Lq N ^ = 1— tjn with tjn — > 0. 
Then, P N -)• if and only if rj N = u(l/VN). 

Proof. For sufficiently large N, 

ViV^ L < log/- 1 < 2^1og ((if')- 1 ) ■ 

We conclude that Pm — > if and only if 

>/tflGg((L<*V)->oo. 

Therefore, 

^Vlog ((ifV 1 ) = -V^log(l - Vn). 

But as a; -» 0, — log(l — x) ~ x/ ln(2). Hence, Pjy — ^ if 
and only if t/nvN —> 00 or rj^ = cj(l/y/N). 

□ 



Now suppose that ci/yN < tj^ < C2/VN. In this case, 
for large iV we deduce that 

ci < logP^ 1 < c 2 , 

or equivalently, 

2" C2 < P/v < 2~ C1 . 

Finally, if 77^ = o(l/y/N) (i.e., r7jv converges to strictly 
faster than l/VN), then P N ->• 1. 

6. CONCLUSION 

We have studied the detection performance of balanced binary 
relay trees with crummy sensors. We have shown that there 
exists an invariant region in the space of (a, f3, q) triplets. We 
have also developed total error probability bounds at the fu- 
sion center as functions of N for both even-height trees and 
odd-height trees. These bounds imply that the total error 
probability converges to sub-linearly, with a decay exponent 
that is essentially y/~N. Compared to balanced binary relay 
trees with no sensor failures, the step-wise shrinkage of the 
total error probability for the crummy sensor case is slower, 
but the asymptotic decay rate is the same. In addition, we 
allow all sensors to be asymptotically bad, in which case we 
deduce necessary and sufficient conditions for the total error 
probability to converge to 0. 
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